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Foreword 


This report documents two new implementations of equivalent linearization for solving geometrically 
nonlinear random vibration problems of complicated structures. The implementations are given the 
acronym ELSTEP, for “Equivalent Linearization using a STiffness Evaluation Procedure.” Both 
implementations of ELSTEP are fundamentally the same in that they use a novel nonlinear stiffness 
evaluation procedure to numerically compute otherwise inaccessible nonlinear stiffness terms from 
commercial finite element programs. The commercial finite element program MSC/NASTRAN 
(NASTRAN) was chosen as the core of ELSTEP. The FORTRAN implementation calculates the 
nonlinear stiffness terms and performs the equivalent linearization analysis outside of NASTRAN. The 
Direct Matrix Abstraction Program (DMAP) implementation performs these operations within 
NASTRAN. Both provide nearly identical results. Within each implementation, two error minimization 
approaches for the equivalent linearization procedure are available — force and strain energy error 
minimization. Sample results for a simply supported rectangular plate are included to illustrate the 
analysis procedure. 

The majority of this work was originally performed by Alexander A. Muravyov during his tenure as a 
National Research Council post-doctoral research associate at the NASA Langley Research Center. 
Additional enhancements and interfaces were subsequently developed by Stephen A. Rizzi. The authors 
wish to thank Jay H. Robinson and Travis L. Turner of the Structural Acoustics Branch at the NASA 
Langley Research Center for helpful discussions and comments. 
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1. Introduction 


Several methods exist for the prediction of geometrically nonlinear dynamic structure response including 
perturbation, Fokker-Plank-Kolmogorov (F-P-K), Monte Carlo simulation and equivalent linearization 
techniques. Perturbation techniques are limited to weak geometric nonlinearities. The F-P-K approach 
[1,2] yields exact solutions, but can only be applied to simple mechanical systems. Monte Carlo 
simulation is the most general method, but computational expense limits its applicability to rather simple 
structures. Finally, equivalent linearization methods [2-6] have seen the broadest application for 
prediction of geometrically nonlinear dynamic response because of their ability to accurately capture the 
response statistics over a wide range of response levels while maintaining a relatively light computational 
burden. 

Implementations of equivalent linearization using finite element analysis have been limited to special 
purpose computer codes. This is largely due to the inaccessibility of the nonlinear stiffness quantities in 
commercial finite element applications. That changed when an equivalent linearization analysis was 
introduced into MSC/NASTRAN (NASTRAN) [7] as a Direct Matrix Abstraction Program (DMAP) 

Alter [8, 9]. In that implementation, the equivalent stiffness was obtained as the sum of the linear 
stiffness and three times the differential stiffness. This formulation was found to over-predict the degree 
of nonlinearity and produce non-conservative results. Over-prediction of nonlinearity can produce the 
undesirable result of structural designs incapable of withstanding the applied loads in an acceptable 
fashion. 

An activity was recently undertaken to more accurately determine the equivalent stiffness through a novel 
approach [10, 1 1], In it, the nonlinear stiffness coefficients from commercial finite element programs 
may be numerically extracted by solving a series of inverse linear and nonlinear static problems. While 
this approach is applicable to any commercial finite element program having a nonlinear analysis 
capability, NASTRAN was selected due to its widespread use in the aerospace industry. The use of this 
new approach in an equivalent linearization analysis has been validated against F-P-K [11] and numerical 
simulation analyses [12] for clamped-clamped beams. 

This report documents two new implementations of the above approach. The implementations are given 
the acronym ELSTEP, for “Equivalent Linearization using a STiffness Evaluation Procedure.” Both 
implementations are fundamentally the same in that they each use the stiffness evaluation procedure 
indicated above. The FORTRAN implementation calculates the nonlinear stiffness terms and performs 
the equivalent linearization analysis outside of NASTRAN. The DMAP implementation performs these 
operations within NASTRAN. Both perform NASTRAN normal modes (solution 103), linear static 
(solution 101) and nonlinear static (solution 106) analyses for the calculation of the nonlinear stiffness 
terms and provide nearly identical results. Optional post-processing is performed by substituting the total 
equivalent stiffness matrix in place of the linear stiffness matrix in a random (solution 111) analysis. 

Within each implementation, two error minimization approaches for the equivalent linearization 
procedure are available - force and strain energy error minimization. Either or both may be run to obtain 
the total equivalent stiffness and root-mean-square displacements. The traditional force error 
minimization approach [3, 4] is implemented as described in [ 10, 12]. An extension of a single degree-of- 
freedom strain energy error minimization approach [5, 6] to multiple degree-of-freedom systems [10] is 
also implemented. 

It should be noted that this analysis has been only developed for, and only validated against cases in 
which the structure exhibits stretching of the middle surface, e.g. clamped and simply supported 
structures. The reader is advised not to apply it to problems in which the nonlinear behavior is manifested 
in other ways, such as in cantilevers. 

Information about the ELSTEP source code distribution and installation may be found in Appendix A. 

The FORTRAN implementation is platform independent. All FORTRAN programs, in both the 
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FORTRAN and DMAP implementations, should be compiled with a FORTRAN 90 compiler. The 
DMAP implementation uses string-based Alters of NASTRAN version 70.0.0 solutions. It is therefore 
expected to be upward compatible with future versions of NASTRAN. 

Throughout this document, specific filenames used in the analysis are made reference to in bold font. 
Sample results for a simply supported rectangular plate are included to illustrate the analysis procedure. 
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2. Finite Element Model Development 


In order to perform the analyses described in Sections 3 - 5, it is necessary to develop three bulk data files 
based on two finite element model databases. One database is for a dynamic model and one is for a static 
model. The dynamic model is used for normal modes (solution 103) and random (solution 111) analyses. 
The static model is a modification of the dynamic model and is used for linear static (solution 101 ) and 
nonlinear static (solution 106) analyses. In the following, it is assumed that the finite element models are 
built in MSC/PATRAN (PATRAN) [13], although this is not required. 


2.1. Normal Modes Analysis Bulk Data File 

• Create a new database called “filename^ g,” for example. Note that in the following, filename is 
an 8-character prefix. 

• Create the finite element model as usual including boundary conditions, material properties, 
loading, etc. 

• In the “Analysis” dialog box, type “filename_e g” as the job name. See for example Figure 1 . 


[\ Analysis jPRjO| 

Action: 

Analyze -j | 

Object: 

Entire Model «j J 

Method 

Analysis Deck -j J 


1 


Code: MS C/ N AST RAN 
Typo: j Struct u ral 


Study:! 


j Available Jobs 

nr 

1 ssptatel.eg j 

| ssptatel.rd 


! fT, 

3 


Job Name 



ppsplatel.og 



Job Description 



MS C/NAS TRAN job created on 
27-Oct-OO at 15:490$ 

p 

j 

tj; 



Translation flirt motor*. .. j j 


Solution Typo... J 

Direct Taxt Input... J 


Subca se Create... | 

Subcase Select... j | 


Apply | 


Figure 1 : Screen capture of PATRAN "Analysis" dialog box for solution 103. 

• In the “Solution Type” dialog box, choose “Normal Modes” analysis (solution 103), as shown in 
Figure 2. 
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Set up the “Solution Parameters’ 
lumped or coupled. 


■ dialog box as shown in Figure 3. The mass calculation may be 


\ Solution Type 


msc/nastran 
S olution Type 


: Solution Typo: 

I LINEAR STATIC 
v NONLINEAR STATIC 
j A NORMAL MOOD 
H v- BUCKLING 
| v COMPLEX EIGENVALUE 
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v TRANSIENT RDPONSE 
j! v NONLINEAR TRANSIENT 
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i 
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OK | 

Cancel | 


Figure 2 


: Screen capture of PATRAN "Solution Type" dialog box for solution 103. 
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r Database Run 



J Cyclic Symmetry 



r Automatic Constraint* 



j Use 5 he II Normals 



Tolerance Anqta 
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1 
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Dynamic Reduction... j 



Figure 3: Screen capture 


°K J Default* | CAtica\ j 

of PATRAN "Solution Parameters" dialog box for solution 103. 
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• From the “Analysis” dialog box, select the “Subcase Create” dialog box and select the 
appropriate load subcase (one with boundary conditions but no loads), as shown in Figure 4. 

• From the “Subcase Parameters” dialog box of the “Subcase Create” dialog box, set the options as 
shown in Figure 5, with the number of desired roots applicable for the particular analysis. This 
set of eigenvectors is subsequently referred to as the “large set of eigenvectors. 



Figure 4: Screen capture of PATRAN "Subcase Create" dialog box for solution 103. 

• From the “Output Requests” dialog box of the “Subcase Create” dialog box, choose the 
“Eigenvectors” result type as shown in Figure 6. 

• Select the “Analysis Deck” method from the “Analysis” dialog box to create the file 
“filename eg. bdf' for subsequent analysis. The extension “_eg” is required as it is referenced in 
the subsequent analyses. An example bulk data file is provided in Appendix B. 
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REAL EIGENVALUE EXTRACTION 


Extraction Method: 


Frequency Range of Interest 


Estimated Number of Rom* 


Number of Desired Roots - 



• Number of Mod»s in Error Analysis 


Figure 5: Screen capture of PATRAN "Subcase Parameters" dialog box for solution 1 03. 


, Output Requests 


SUBCASE NAME: Default 

SOLUTION SEQUENCE: 103 




Select Result Type 
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; Output Requests 
! ( J VECTORffORTI ,REAL)=AH FEM 


Defaults Cancel 


Figure 6: Screen capture of PATRAN "Output Requests" dialog box for solution 103. 
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2.2. Static Analysis Bulk Data File 

• Copy the dynamic database “filename _eg.db" to “filename jt.dbf for example. The 8-character 
prefix filename must be the same as that used for the normal modes bulk data file. 

• Open “filename _st. db“ in PATRAN and make the following modifications: 

o Add an extra grid point with a connecting beam or rod element. Note the added element 
must have a nonlinear capability and be connected to the structure at a fixed boundary. 
Equivalence the node where the new element is connected to the structure and, if 
necessary, renumber so that the node numbers of the main structure are the same as those 
for the normal modes analysis. In this manner, the extra grid point and element number 
will be greater than the highest grid point and element number of the original model. 

o Create a new load case and delete any previously applied loads from the model. To the 
new load case, apply an arbitrary non-zero nodal force at the extra grid point. 

o To the new load case, create a constraint displacement set, which sets all degrees of 
freedom to zero values except for those at the extra grid point. 


Note that the necessity to introduce an extra element is an artifact of the solver 
implemented in NASTRAN solutions 101 (linear static) and 106 (nonlinear static). 
Specifically, these solutions can only solve the forward static problem, that is, these 
solutions solve for a set of displacements from a specified loading. In order to determine 
the nonlinear stiffness coefficients, it is necessary to solve linear and nonlinear inverse 
problems, which compute the forces due to a prescribed set of displacements. The extra 
element is introduced to allow NASTRAN to solve for a dummy set ot displacements 
with the set of prescribed displacements acting as displacement constraints. The sought 
vector of applied loads (not including the extra load) is identical to the set of single-point- 
constraint forces in NASTRAN terminology. 


\ Solution Parameters 


j Nonlinear Static Solution Parameter* 


J Automatic Constraints 
r Urge Displacements 
J follower Forces 
Mass Calculation: 

Data Deck Echo: 

Plate Rat Stiffness Factor = 
Maximum Printed Lines = 
Maximum Run Time = 
Wt.-Mass Conversion = 
Node i d. for Wt. Goner. = 


Coupled -i | 
None -j j 


1999999999 


prr 


OK 


Defaults [ Cancel 


Figure 7: Screen capture of PATRAN "Solution Parameters" dialog box for solution 106. 


• In the “Analysis” dialog box, type “filename jX" as the job name. 
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• In the “Solution Type” dialog box, choose “Nonlinear Static” analysis (solution 106). Note that 
the resulting “.bdf” file is also used in solution 101 and is automatically modified for this 
purpose, as described in Sections 3 and 4. 

• Set up the “Solution Parameters” dialog box as shown in Figure 7. The mass calculation is 
irrelevant for the static case. 

• From the “Analysis” dialog box, select the “Subcase Create” dialog box and select the new load 
subcase. 

• From the “Subcase Parameters” dialog box of the “Subcase Create” dialog box, set the options as 
shown in Figure 8. 



Figure 8: Screen capture of PATRAN "Subcase Parameters" dialog box for solution 106. 

• From the “Output Requests” dialog box of the “Subcase Create” dialog box, choose the 
“Constraint Forces” result type as shown in Figure 9. 

• Select the “Analysis Deck” method from the “Analysis” dialog box to create the file 
"filename _st. bdf for subsequent analysis. The extension “_st” is required as it is referenced in 
the subsequent analyses. An example bulk data file is provided in Appendix C. 

2.3. Random Analysis Bulk Data File 

A random analysis (solution 1 1 1) is required only if post-processing of displacement power spectral 
densities or stress/strain root-mean-square values or PSDs is desired. 

• In the “ filename _eg.db" database, create a new load subcase and add the appropriate dynamic 
loading (in physical coordinates) corresponding to the modal loading specified in the equivalent 
linearization analysis (see Sections 3. 1 .2 and 4. 1 .2). 
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Figure 10: Screen capture of PATRAN "Solution Parameters" dialog box for solution 111. 
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• In the “Solution Type” dialog box, choose “Frequency Response” analysis. 

• In the “Solution Parameters” dialog box, select the settings as indicated in Figure 10 with the 
formulation as modal to specify solution 1 1 1 and mass calculation as desired (lumped or 
coupled). 

• Select the “Eigenvalue Extraction” dialog box from the “Solution Parameters” dialog box and set 
up the values as indicated in Figure 1 1. The number of desired roots should be the same as the 
large set of eigenvectors specified in the development of the normal modes analysis bulk data file. 



REAL EIGENVALUE EXTRACTION 



Extraction Method: La not os ~J j 

Frequency Range of Interest 


| Lower = p 

] Upper = p. 


Estimated Number of Roots 


pi 00 

Number of Desired Roots = 



Figure 1 1 : Screen capture of PATRAN "Eigenvalue Extraction" dialog box for solution 111. 

• From the “Subcase Create” dialog box, specify the frequency range in the “Subcase Parameters” 
dialog box. See Figure 12, for example. For consistent results between the equivalent 
linearization analysis and the solution 1 1 1 post-processing analysis, the frequency range should 
be the same as that used in the equivalent linearization analysis (see Sections 3.1.4 and 4. 1 .4). 


[\ Subcase Parameters 



Starting Frequency = 

pp.o 


Ending Frequency = 

jl02< 


# of Freq. Increments = 

1 4096 

1 ! 

°* 1 


Cancel [ 


Figure 12: Screen capture of PATRAN "Subcase Parameters" dialog box for solution 111. 

• In the “Analysis” dialog box, type fixeDlll_rd as the job name to create the bulk data file 
fixeDl 1 l_rd.bdf. This filename is required as it is referenced in the subsequent analysis. 

• It is necessary to manually edit the bulk data file fixeDl 1l_rd.bdf in preparation for the analysis 
run. The manual edits are indicated in bold in the sample listing provided in Appendix D. 
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3. FORTRAN Implementation 


The analysis is subdivided into two parts. The first part calculates the nonlinear stiffness coefficient^ 

The second part performs one of two equivalent linearization procedures to compute the total equi 
stiffness matrix and RMS displacements. A description of optional post-processing is provided in Section 

5. 


3.1. Input Files 

In addition to the bulk data files discussed in Section 2, several additional input files are required to 
specify various parameters used in the analysis. These are detailed below. 

3.1.1. Modes Selection File 


The file fixed.mod indicates which eigenvectors out of the large set of eigenvectors participate m the 
analysis. The format for fixed.mod is a free ASCII format as specified below. 

LINE 1 : Number of selected modes (NMOD). 


LINE 2: Selected modes (e.g. I s ' and 4' h ). 

LINE 3: Scaling coefficients for each eigenvector selected. The product of each eigenvector and scaling 
coefficient produces a displacement field, which is used as a prescribed displacement set. he 
scaling coefficients should be chosen to be equal and not made too large so that the product of 
the eigenvector and coefficient will represent a realistic prescribed displacement within t e 
linear range. The value of 1 .Oe-4 appears to work well for cases considered thus far. 


A sample listing of fixed.mod with two selected modes (NMOD-2) is provided below. 


2 


1 


4 


1 . e- 04 1 . e-04 

Note that specification of the participating modes may necessitate running a standard normal modes 
analysis prior to the procedure for calculating the nonlinear stiffness coefficients described in Section 3.2. 
In addition, if a post-processing analysis is to be performed (as described in Section 5), the same 
information must also be specified in the file fixeD.mod, as described in Section 4.1.1. 

3.1.2. Spectral Density Loads File 


The file fixed.den contains the spectral density matrix of the applied loads in modal coordinates ~S a (co) . 
The format for fixed.den is a free ASCII format as specified below. 

LINE 1 : Number of breakpoints used to define frequency spectrum. (Min 2, Max 50). Between each 
breakpoint, each value of the Sjr (to) is linearly interpolated at the frequency increment 
specified in fixed.par (see Section 3.1 .4). 

LINE 2: Frequency of first breakpoint (Hz) 

LINE 3: First row of modal spectral density matnx (NMOD entries) at first frequency breakpoint. 


LINE Last row of modal spectral density matrix (NMOD entries) at first frequency breakpoint 
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(Lines 2 to Line 3+NMOD are repeated for each breakpoint) 


A sample listing of fixed.den for a flat spectrum between 0-1024 Hz with two selected modes is provided 
below. 


2 

0 . 00000000000000 
6 . 83899171960566 
2 .27952175755446 
1024 . 00000000000 
6 . 83899171960566 
2 . 27952175755446 


2.27952175755947 
0 . 759793205814123 

2 .27952175755947 
0 . 759793205814123 


To be explicit, at both 0 and 1024 Hz, the following spectral density matrix of applied loads in modal 
coordinates is specified: 


Sjr (0,1024) 


Sffn 

Sffn 


'6.84 

2.28' 

_Sff 21 

Sff 22 _ 


2.28 

0.76 


This matrix is computed as 

Sff{ai) = <t> T S ir {a))<l) 0) 

where S g {(o) is the fully-populated, real-valued spectral density matrix of the load in physical 
coordinates and <p is the subset of selected eigenvectors. Because S ff (ft>) is real-valued, only correlated 

inputs are permitted. This condition allows loadings such as normal acoustic pressure, inertial loading 
(for base excitation), and one or more correlated point forces or moments. It does not allow loadings 
which generate complex terms, such as grazing incidence acoustic loading. 


The values of S g (o>) specify the double-sided spectrum level in units7rad/sec. In order to specify the 
same load in physical coordinates for a solution 1 1 1 post-processing analysis, the values of S ff {(£>) must 

be converted to a single-sided spectral density in units 2 /Hz by multiplying by 2 x 2 tc. These may then be 
specified on the TABRND1 card of the bulk data file. 

From equation (1), it is obvious that a normal modes analysis must first be performed to obtain the 
eigenvectors. Several files containing eigenvectors are formed in the first part of the analysis used to 
calculate the nonlinear stiffness coefficients (see Section 3.2). The eigenvectors to be used in equation 
( 1 ) are the set of selected g-size (all degrees of freedom including constrained ones) eigenvectors from the 
file egveCa.dat. 

In order to facilitate generation of the modal spectral density matrix of the applied loading, the utility 
program gensff (described in Section 3.3) may be used to automatically compute the values and write 
them to fixed.den for a limited set of loading conditions. 
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3.1.3. Damping Matrix File 

The file fixed.dam contains the diagonal modal damping matrix. The format for fixed.dam is a free 
ASCII formatted square matrix (NMOD x NMOD) as specified below. 


c, o ••• o 

0 C 2 0 0 

0 C NMOD 


A sample listing of fixed.dam with two selected modes is provided below. 

14.6618 0 . 

0 . 14.6618 

In order to specify the same damping for a solution 1 1 1 post-processing analysis, the values must be 
converted to percent critical damping through the usual relation 


where C are the modal damping values (as specified in fixed.dam) and (O, are the linear eigenvalues, not 
those of the equivalent linear system. This value of load may then be specified on the TABDMP1 card of 
the bulk data file. 

3.L4. Parameter File 

The file fixed.par contains various additional parameters required for the analysis. The format tor 
fixed.par is a free ASCII format as specified below. 


LINE 1 : Minimum and maximum frequency range of the analysis (Hz). Note that the ™ xl1 ™ 

frequency range should be selected well beyond the frequency of the highest selected mode. 

This is necessary so that the resonant frequencies don’t shift out of the analysis bandwidth as 
nonlinearity increases. If the minimum or maximum analysis frequencies specified above fall 
outside of the bandwidth of applied loads (specified in fixed.den), the load will be padded with 
zeros at either low or high frequency end. If the minimum or maximum analysis frequencies fall 
inside of the bandw idth of applied loads, the loading will start and end at the frequencies 

specified above. 


LINE 2: Frequency increment used in the analysis (Hz) 

LINE 3: Weightings a and p . A discussion on how the weightings are used is given in [12]. The 
weightings must sum to 1 . Typical initial values are Of = 0.4 and = 0.6 for the force error 
minimization approach and a = 0.05 and (3 = 0.95 for the strain energy error minimization 
approach. The choice of weightings may influence the solution to which the system converges, 
as discussed in Section 3.3.1. 


LINE 4: Convergence criteria £ 


A sample listing of fixed.par is provided below. 
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0 . 1024 . 

0 .25 

0.4 0.6 

0 .005 


3.2. Calculation of Nonlinear Stiffness Coefficients 

The method for calculating nonlinear stiffness coefficients is shown in the flowchart in Figure 1 3. Each 
operation is separated by a dashed line and is run in the sequence shown from top to bottom. The vano 
operations are described below. 

• nreDl03 - This FORTRAN program reads the bulk data file “ filename _eg.bdf and rewrites it to 
the file fixed 103_eg.bdf. The program prepl03 adds lines to the file management and case 
control sections of the bulk data file ‘filename _eg.bdf" The lines that are added are shown in 
bold in Appendix E. It also creates the file prob_prefix.scr, which contains the 8-character 
prefix (i.e. filename) for use in other file manipulation programs. 

• The first NASTRAN run calculates the eigenvectors of the model and writes them to two Files for 
subsequent use as described below. This NASTRAN run utilizes the DMAP alter dm_103.v to 
modify the standard normal modes analysis (solution 103). 

o egvec.dat contains the mass-normalized large set of g-size eigenvectors (all degrees of 
freedom including constrained ones) requested in the modeling phase (see Figure - ). 

This is an ASCII file with the INPUTT4 DMAP matrix input format [14]. 

o egveH.dat contains the mass-normalized large set of 1-size eigenvectors (all 

unconstrained degrees of freedom) requested in the modeling phase. This is a binary file. 

• gendisp - This FORTRAN program generates the displacements fields used in the subsequent 
solution 106 and 101 runs. The following files are created: 

o egveCa.dat contains the mass-normalized selected set of g-size eigenvectors specified in 
the file fixed.mod. This data is a subset of egvec.dat. This is an ASCII file with the 
INPUTT4 DMAP matrix input format. 

o displ.pr2 contains the number of unique combinations of two modes (IN) and the 

combinations of those modes used to prescribe the displacement fields imposed for the 
inverse problems. This is an ASCII file in free format. 

o displ.pr3 contains the number of unique combinations of three modes (IC) and the 
combinations of those modes used to prescribe the displacement fields imposed for the 
inverse problems. This is an ASCII file in free format. 

o dispL inn are the prescribed displacement fields used for the linear inverse problem. 
This is an ASCII file with the INPUTT4 DMAP matrix input format. 

o displ.inn are the prescribed displacement fields used for the nonlinear inverse problem. 
This is an ASCII file with the INPUTT4 DMAP matrix input format. 

• Drepl06 - This FORTRAN program reads the bulk data file ‘filename _st.bdf and rewrites it to 
the file fixedl06_st.bdf. The program prepl06 adds lines to the file management and case 
control sections of the bulk data file “filename jt.bdf-' In particular, it writes the number osu 
cases (i e. the number of prescribed displacement fields) required in the subsequent NASTRAN 
run The number of sub-cases depends on the number of modes selected (NMOD) number of 
unique combinations of two modes (IN) from the number of modes selected, and the number 
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unique combinations of three modes (IC) from the number of modes selected. It is calculated 
using the following relation 

NSUBC = 2xNMOD+3xlN + IC 






Thus for NMOD = 3, IN = 3 (combinations of modes 1&2, 2&3, and 1&3), IC = 1 (combination 
of modes 1&2&3) and the number of sub-cases NSUBC =16. As an example, a listing of the 
changed lines of file fixed!06_st.bdf are shown in bold in Appendix F for NMOD - 2. 

The second NASTRAN run calculates the force vector corresponding to the prescribed 
displacement fields (provided in the file displ.inn) in a nonlinear analysis. Because it is not 
possible to separate the linear and nonlinear components, the resulting force vector contains the 
combined force. The resulting force vector is written to f_N.frc, an ASCII file with the 
INPUTT4 DMAP matrix input format. The DMAP alter dm_106.v is used to modify the 
standard nonlinear static analysis (solution 106). 

preplOl - This FORTRAN program reads the bulk data file “filename s/. bdf and rewrites it to 
the file FixedlOlst.bdf. The program preplOl adds lines to the file management and case 
control sections of the bulk data file “ filename _st. bdf. ” In particular, it writes the number of sub 
cases (i e the number of prescribed displacement fields) required in the subsequent NASTRAN 
run The number of sub-cases is equal to the number of modes selected (NMOD). As an 
example, a listing of the changed lines of file fixed 101_st.bdf are shown in bold in Appendix G 

for NMOD = 2. 


The third NASTRAN run calculates the force vector corresponding to the prescribed 
displacement fields (provided in the file dispL.inn) in linear analysis. The resulting force vector 
is written to flL.frc, an ASCII file with the INPUTT4 DMAP matrix input format. The DMAP 
alter dm 101 . v is used to modify the standard linear static analysis (solution 101). 


nlcfs - This FORTRAN program calculates the nonlinear stiffness coefficients of the structure 
and stores them in the file nlcfS.dat. Note that both quadratic a jk and cubic b jU stiffness 

coefficients [10-12] are computed in nlcfs, but only the cubic terms are stored in nlcfS.dat 
because a zero mean response is assumed. The file nlcfS.dat is an ASCII file with the 
FORTRAN format (1P,3E23.16). 


The nonlinear stiffness coefficients computed through this series of operations depend only on the 
structure (geometry, material properties, and boundary conditions) and not on the loading level or 
distribution. Therefore, this portion of the analysis only needs to be performed once for a given structural 
configuration. The effect of loading is addressed in the next sequence of analyses (see Section 3.3). 

In order to check the functioning of the entire analysis procedure, it is recommended that a linear analysis 
be run from this point forward using the two remaining sequences of analyses (see Sections 3.3 and 5 as 
is. To do so, all nonlinear stiffness coefficients in nlcfS.dat should be set to very small (non-zero) values 
(e.g. 10' 15 ), and the weighting a and /? in the file fixed.par should be set to 1 and 0, respectively. If 
everything is working well, the results obtained will be virtually identical to those from a linear solution 
1 1 1 analysis. 


The UNIX script file doitl.bat (see Appendix A) runs the first part of the analysis as outlined above in an 
automated fashion. Some commands in doitl.bat may need to be modified depending on the operating 
system and how NASTRAN is installed. Removal (or transfer to another directory) of previously created 
ELSTEP output files is necessary for proper program execution. Removal of previously created standard 
NASTRAN files, e g. *.fl)6, *.fD4, *.op2, *.log, and *.pch, is optional, but helps to reduced clutter. 
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Figure 13: Flow-chart for calculation of nonlinear stiffness coefficients in the FORTRAN 

implementation. 
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3 . 3 . Equivalent Linearization Procedure 


shown from top to bottom. The vanous operat.ons are descnbed below. 


OPTIONAL 



Figure 14: 


Flow chart for calculation of total equivalent stiffness matrix and RMS displacements m the 

FORTRAN implementation 




ppnsff-The first operation is optional and is used to develop the frequency dependent spectral 

specify an inertial load) and point forces and moments- In each case, the loadings are u y 

"heZ^ng is read from fixed..!!. The computed modal spectral denstry matnx ,s wntten to 
fixed.den (as described in Section 3.1.2). 

o The parameter file fixed.sff specifies the following information in free format: 


LINE 1 : Load option (LO ADOPT) . , » 

1 = uniformly distributed pressure with uniform mesh 

2 = point force(s) and moment(s) 


If LO ADOPT =1 

LINE 2: AREA, NOD, NODB 

AREA - Surface area of structure 

NOD - Total number of nodes in the model 
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NODB - Number of constrained (boundary) nodes 

Forces at the grid points are obtained by multiplying the pressure level by the 
area and distributing these over the mesh, hence the need for a uniform mesh. If 
a non-uniform mesh is used, a pressure loading may be specified using 
LOADOPT = 2 with forces input at each grid point. 

LINE 3: NPT 

NPT _ Number of breakpoints in the PSD ( 2 < NPT < 50 ) 

LINE 4: (repeated NPT times): FREQ, PSD 

FREQ, PSD - Breakpoint pairs for PSD definition of uniformly distributed 
pressure. Values of PSD will be linearly interpolated between breakpoint pairs 
at the frequency increment specified in fixed.par. FREQ is the frequency 
specified in Hz. PSD is the single-sided pressure PSD in units /Hz and is the 
same as would be used on the TABRND1 card for a solution 1 1 1 analysis. Note 
that gensff converts the load internally to a double-sided value in units /rad/sec. 

Note that in this case of a uniformly distributed pressure on a uniformly 
distributed mesh, the forces for each grid point are identical and act in the z- 
direction. Hence, all non-zero auto- and cross-spectral density values in S ff are 
equal. Further, each column is identical. These properties speed the 
computation of S ff - . 

If LOADOPT = 2 
LINE 2: NLOAD 

NLOAD Number of fully correlated point forces and moments. 

LINE 3: NPT 

NPT - Number of breakpoints in the PSD ( 2 < NPT < 50 ) 

(Lines 4-5 repeated NLOAD times) 

LINE 4: NODF, DOF1-DOF6 

NODF - Grid point at which point force is applied 

DOF 1 -DOF6 — Multiplier for each of six degree-of-freedoms at each grid point. 
An arbitrarily oriented force or moment may be applied in this fashion. 

LINE 5: (repeated NPT times): FREQ, PSDF 

FREQ, PSDF - Breakpoint pairs for PSD definition of point forces. Values of 
PSDF will be linearly interpolated between breakpoint pairs at the frequency 
increment specified in fixed.par. FREQ is the frequency in Hz. PSDF is the 
single-sided force PSD in units 2 /Hz and is the same as would be used on the 
TABRND1 card for a solution 1 1 1 analysis. Note that gensff converts the load 
internally to a double-sided value in units 2 /rad/sec. 

Note that in this case, the forces for each grid point are not identical. Therefore, 
the non-zero auto- and cross-spectral density values in S ff are unequal and each 

column is not identical. These properties make the computation of S ff more 

intensive. 
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Note that in general, the cross-spectral densities of the input loading are required to 
specify S g . However, since both load options require fully correlated (in phase) 
loadings, the cross-spectral density (CSD) terms may be obtained from the auto-spectral 
density (PSD) from the following relationship: 

CSD,, = ±^{PSD,)(PSD) 

This allows for the specification of only the PSD values, and not both the PSD and CSD 
values, as indicated above. Note that gensff as currently implemented only takes the 
positive square root. Therefore, the specified loadings must have the same direction. 
Failure to do so will result in an inconsistency in the sign of the diagonal and off-diagonal 

terms in S g - . 

In the upcoming examples, the following may be helpful. The values specified above for the 
spectral density matrix of the loading in physical coordinates ( S lf ) may be determined from an 
RMS level for a particular bandwidth using the following relation 

RMS 2 / 

S = — — 

JJ 2x2n 

where B W is the frequency bandwidth in Hz. For example, for an RMS acoustic pressure of 2048 
Pa (160.21 dB) and 0-1024 Hz bandwidth, the double-sided value of S ff is 325.95 Pa /rad/sec. 
The corresponding single-sided value, which should be used in fixed.sff, is 4096 Pa /Flz. 


The following is an example of fixed.sff for the uniformly distributed pressure indicated above: 
l 

0.0903224 609 96 
0 . 4096 . 

1024 . 4096 . 


The following is an example of fixed.sff for point forces. The force at grid point 305 acts in the 
z-direction and linearly increases in magnitude from 0 to 5 N 2 /Hz over the bandwidth 0-512 Hz. 

A second force acts in the x-direction at grid point 220 with a constant magnitude of 2 N'/Hz over 
the frequency range of 0-5 1 2 Hz. 

2 
2 
2 

305 0.0 

0 . 0.0 

512 . 5.0 

220 2.0 

0 . 1.0 


0.0 1.0 0.0 0.0 0.0 


0.0 0.0 0.0 0.0 0.0 
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512 . 


1 . 0 

Note that when converting from a uniform pressure load specified by load option 1 to a force load 
specified by load option 2, the PSD values specified in option 1 must be multiplied by the area 

squared. 

• prepel - This FORTRAN program reads the NASTRAN iD6 file fixedl03_eg.f06 produced in 
the second step of the previous procedure (see Section 3.2) and writes the selected natural 
frequencies to the file fixed.stf. This is an ASCII file in free format. 

• The third operation performs the equivalent linearization procedure through the FORTRAN 
programs eqlf and eqls, for the force error and strain energy error minimization approaches, 
respectively. Each of these alters produce the following output. 

o st_f.dat and st_s.dat contain the total equivalent stiffness matrix in modal coordinates for 
the force and strain energy error minimization approaches, respectively. These are ASCII 
files with the INPUTT4 DMAP matrix input format. 

o rms_f.dat and rms_s.dat contain the root-mean-square displacements in physical 
coordinates for all degrees of freedom. These are ASCII files in free format. 

o dcov f.dat and dcov_s.dat contain the cross-covariance matrix of modal displacements. 
These are ASCII files in free format. 

The UNIX script files doit2_f.bat and doit2_s.bat (see Appendix A) run the second part of the analysis 
as outlined above in an automated fashion. 

3.3.1. Convergence 

The weightings a and [3 . and the convergence criteria £ , play a significant role in the convergence of 
the equivalent linearization procedure. Because the system is nonlinear, it is possible for multiple 
solutions to exist for a given loading condition. The selection of different weightings may cause the 
equivalent linearization procedure to converge to different solutions. Therefore, it is advisable to explore 
the solution space by trying out different weightings before selection of a converged solution. If the root- 
mean-square displacements from the equivalent linearization procedure, as reported in rms_f.dat or 
rms_s.dat, differ with the choice of weighting, then further analyses may be warranted, such as 
observation of the deflected shape or the equivalent linear power spectral density from a post-processing 
analysis (see Section 5). 
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4. DMAP Implementation 


As in the FORTRAN implementation, the analysis for the DMAP implementation is subdivided into two 
parts The first part calculates the nonlinear stiffness coefficients. The second part performs one of two 
equivalent linearization procedures to compute the total equivalent stiffness matrix and RMS 
displacements. A description of optional post-processing is provided in Section 5. 


Note that since some 
implementation, both 


files bear the same name but have a different format from those in the FORTRAN 
FORTRAN and DMAP implementations should not be run from the same directory. 


4.1. Input Files 

In addition to the bulk data files discussed in Section 2, several additional input files are required to 
specify various parameters used in the analysis. These are detailed below. 

4.1.1. Modes Selection File 


As in the FORTRAN implementation, the file fixeD.mod indicates which eigenvectors out of the large set 
of eigenvecion, partiup.« in .he analysis. The forma, for fixeD.mod is a INPUTT4 DMAP mafnx mpo. 
format as specified below. 

LINE 1 • Number of columns (always 2), number of selected modes (NMOD), matrix form (always 2 for 
rectangular), type of matrix (always 2 for real, double precision), name of the matrix (always 
PHG), FORTRAN format specification (always 1P,3E23.16). 

LINE 2: Column number (always 1), row position of first nonzero term (always 1), number of real 
double-precision entries on the following line (NMOD). 

LINE 3: Selected modes (up to three per line). Note that these numbers are input as real and converted to 


integers within the program 

LINE 4: Column number (always 2), row position of first nonzero term (always 1), number of real 
double-precision entries on the following line (NMOD). 

LINE 5: Scaling coefficients for each eigenvector selected (up to three per line). See Section 3. 1 . 1 for a 
discussion of these terms. 


A sample listing of fixeD.mod with two selected modes (NMOD-2) is provided below. 
2 22 2PHG IP , 3E23 . 16 

1 1 2 

1 . OOOOOOOOOOOOOOOOE+OO 4. OOOOOOOOOOOOOOOOE + OO 
2 1 2 

1. 000000000000000 0E- 04 1 . 000000000000000 0E- 04 

It is clear that the size of this file depends on the number of modes selected. Since the fonmt ain be 
somewhat confusing, a sample listing with four selected modes (NMOD=4) is also provided below. 


2 PHG 


IP, 3E23 . 16 


2 4 2 

1 1 4 

1.0000000000000000E+00 3 . OOOOOOOOOOOOOOOOE+OO 7 . OOOOOOOOOOOOOOOOE+OO 
1 . 0000000 OOOOOOOOOE+Ol 
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2 


1 


4 


1 .OOOOOOOOOOOOOOOOE-04 1 . 0000000000 OOOOOOE- 04 1 . OOOOOOOOOOOOOOOOE - 04 
1.0000000000000000E-04 

4. 1.2. Spectral Density Loads File 

The file fixeD.den is analogous to fixed.den in the FORTRAN implementation, except that the spectral 
density matrix of the applied loads in modal coordinates S j is independent of frequency in the DMAP 
implementation, i.e. only flat spectra may be applied. The format for fixeD.den is a INPUTT4 DMAP 
matrix input format as specified below. 

LINE 1 : Number of selected modes (NMOD), number of selected modes (NMOD), matrix form (always 

1 for square), type of matrix (always 4 for complex, double precision), name of the matrix 
(always PHG), FORTRAN format specification (always 1P,6E13.6). 

LINE 2: Column number (always 1), row position of first nonzero term (always 1), number of real 
double-precision entries on the following line (2 x NMOD). 

LINE 3: First column of modal spectral density matrix (NMOD entries, up to 3 complex entries per line). 
Each entry is complex with the imaginary part equal to zero. 

(Lines 2 and 3 repeat for NMOD selected modes) 

2 nd TO LAST LINE: 

Column number (always NMOD), row position of first nonzero term (always 1), number of real 
double-precision entries on the following line (2 x NMOD). 

LAST LINE: 

Last column of modal spectral density matrix (NMOD entries, up to 3 complex entries per line). 

A sample listing of fixeD.den for two selected modes is provided below. 

221 4 PHG IP, 6E13 . 6 

114 

0 . 683899E+01 O.OOOOOOE+OO 0.227952E+01 0.000000E+00 

2 14 

0 . 227 952E+ 0 1 O.OOOOOOE + OO 0.759793E + 00 0.000000E+00 

To be explicit, the above specifies the following frequency-independent spectral density matrix of applied 
loads in modal coordinates: 


1 

:? 

Coi 

3 

to 

1 


‘(6.84,0) (2.28,0)' 

_*SV 21 S ff 22 _ 


(2.28,0) (0.76,0) 


This matrix is computed as 

~S ff = </> T SJ (2) 

where S ff is the real-valued spectral density matrix of the load in physical coordinates (double-sided 
spectrum level in units 2 /rad/sec) and <p is the subset of selected eigenvectors. The eigenvectors to be used 
in equation (2) are the set of selected g-size (all degrees of freedom including constrained ones) 
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eigenvectors from the file egveCa.dat (see Section 4.2). As in the FORTRAN implementation, only 
correlated inputs are permitted. Values of S must be converted to single-sided values as indicated in 

Section 3.1.2 for use in a solution 1 1 1 post-processing analysis. 

In order to facilitate generation of the modal spectral density matrix, the utility program gensffD 
(described in Section 4.3) may be used to automatically compute the values and write them to fixeD.den 
for a limited set of loading conditions. 

4. 7.3. Damping Matrix File 

The file fixeD.dam is analogous to fixed. dam in the FORTRAN implementation. It contains the 
diagonal modal damping matrix. The format for fixeD.dam is a INPUTT4 DMAP matrix input format as 
specified below. 

LINE 1 : Number of selected modes (NMOD), number of selected modes (NMOD), matrix form (always 

1 for square), type of matrix (always 2 for real, double precision), name of the matrix (always 
PHG), FORTRAN format specification (always 1P,3E23.16). 

LINE 2: Row and column position of first modal damping values (always 1,1), number of real double- 
precision entries on the following line (always 1). 

LINE 3: First modal damping value 

(Lines 2 and 3 repeat for NMOD selected modes) 

2 nd TO LAST LINE: 

Row and column position of last modal damping values (always NMOD, NMOD), number of 
real double-precision entries on the following line (always 1). 

LAST LINE: 

Last modal damping value. 

A sample listing of fixeD.dam with two selected modes is provided below. 

2 2 1 2 PHG 1 P , 3E23 . 16 

111 

0 . 14 6 6 18 0000 000000E+02 
2 2 1 

0 . 1466180000000000E + 02 

In order to specify the same damping for a solution 1 1 1 post-processing analysis, the values must be 
converted to percent critical damping as discussed in Section 3.1.3. 

4.1.4 . Parameter File 

The file fixeD.par is analogous to fixed.par in the FORTRAN implementation. It contains various 
additional parameters required for the analysis. The format for fixeD.par is a INPUTT4 DMAP matrix 
input format as specified below. 

LINE 1 : Number of columns (always 1), number of rows (always 6), matrix form (always 2 for 

rectangular), type of matrix (always 2 for real, double precision), name of the matrix (always 
PHG), FORTRAN format specification (always 1P,3E23.16). 

LINE 2: Column number (always 1), row position of first nonzero term (always 1 ), number of real 
double-precision entries on the following line (always 6). 
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LINE 3: Minimum and maximum frequency range of the analysis (Hz), frequency increment used in the 
analysis (Hz). Note that the maximum frequency range should be selected well beyond the 
frequency of the highest selected mode. This is necessary so that the resonant frequencies don’t 
shift out of the analysis bandwidth as nonlinearity increases. 

LrNE 4: Weightings cc and /? , and convergence criteria £ (see Section 3.1.4 for description). Note that 
these values are put on a fourth line because the format on line 1 indicates up to three values per 
line. The choice of weightings may influence the solution to which the system converges, as 
discussed in Section 3.3.1 . 

A sample listing of fixeD.par is provided below. 

162 2PHG IP, 3E23 . 16 

116 

O.OOOOOOOOOOOOOOOOE+OO 1.0240000000000000E+03 0.2500000000000000E+00 
0.4000000000000000E+00 0.6000000000000000E+00 0.5000000000000000E-02 

4.2. Calculation of Nonlinear Stiffness Coefficients 

The method of calculating the nonlinear stiffness coefficients is shown in the flowchart in Figure 15. 

Each operation is separated by a dashed line and is run in the sequence shown from top to bottom. The 

various operations are described below. 

• prepHBD - This FORTRAN program is the analogue to prepl03 in the FORTRAN 
implementation. It reads the bulk data file “filename _eg.bdf' and rewrites it to the file 
fixeD103_eg.bdf. The program prepHBD adds lines to the file management and case control 
sections of the bulk data file “filename eg.bdf." The lines that are added are shown in bold in 
Appendix H. It also creates the file prob_prefix.scr, which contains the 8-character prefix (i.e. 
filename) for use in other file manipulation programs. 

• The first NASTRAN run is analogous to the FORTRAN implementation, except that the DMAP 
alter dm_103D.v also generates the displacement fields. These were generated by gendisp in the 
FORTRAN implementation. The following files are created: 

o egvec.dat contains the mass-normalized large set of g-size eigenvectors (all degrees of 
freedom including constrained ones) requested in the modeling phase (see Figure 5). It 
contains the same information as that produced in the FORTRAN implementation, but in 
a binary format. 

o egveH.dat contains the mass-normalized large set of 1-size eigenvectors (all 

unconstrained degrees of freedom) requested in the modeling phase. This binary file is 
identical to that produced in the FORTRAN implementation. 

o egveCa.dat contains the mass-normalized selected set of g-size eigenvectors specified in 
the file fixeD.mod. It is identical to that produced in the FORTRAN implementation. 

o egveC.dat is a binary version of egveCa.dat. There is no analogue in the FORTRAN 
implementation. 

o displ.pr2 contains the number of unique combinations of two modes (IN) and the 

combinations of those modes used to prescribe the displacement fields imposed for the 
inverse problems. It contains the same information as that produced in the FORTRAN 
implementation, but in the INPUTT4 DMAP matrix input format. 
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o displ.pr3 contains the number of unique combinations of three modes (IC) and the 
combinations of those modes used to prescribe the displacement fields imposed for the 
inverse problems. It contains the same information as that produced in the FORTRAN 
implementation, but in the INPUTT4 DMAP matrix input format. 

o displ.pr4 contains additional information used to prescribe the displacement fields 
imposed for the inverse problems. This is an ASCII file with the INPUTT4 DMAP 
matrix input format. There is no analogue in the FORTRAN implementation. 

o dispL.inn are the prescribed displacement fields used for the linear inverse problem. It 
contains the same information as that produced in the FORTRAN implementation, but in 
a binary format. 

o displ.inn are the prescribed displacement fields used for the nonlinear inverse problem. 

It contains the same information as that produced in the FORTRAN implementation, but 
in a binary format. 

• prep!06D - This FORTRAN program is the analogue to prepl06 in the FORTRAN 
implementation. It reads the bulk data file “filename _st.bdf' and rewrites it to the file 
fixeD106_st.bdf. The program prepl06D adds lines to the file management and case control 
sections of the bulk data file “filename _st.bdf” As an example, a listing of the changed lines of 
file fixeD106_st.bdf are shown in bold in Appendix I for NMOD = 2. 

• The second NASTRAN run is analogous to the FORTRAN implementation. In the DMAP 
implementation, the file fJV.frc is binary. The DMAP alter dm_106D.v is used to modify the 
standard nonlinear static analysis (solution 106). 

• preplOlD - This FORTRAN program is the analogue to preplOl in the FORTRAN 
implementation. It reads the bulk data file “filename _st.bdf ? and rewrites it to the file 
fixeD101_st.bdf. The program preplOlD adds lines to the file management and case control 
sections of the bulk data file “filename jst.hdf” As an example, a listing of the changed lines of 
file fixeDlOlst.bdf are shown in bold in Appendix J for NMOD = 2. 

• The third NASTRAN run is analogous to the FORTRAN implementation. In the DMAP 
implementation, the file fl_L.frc is binary. The DMAP alter dmlOlD.v is used to modify the 
standard linear static analysis (solution 101). 

• The fourth NASTRAN run is the DMAP analogue of the FORTRAN program nkfs. The file 
nlcfs.bdf is a standalone DMAP alter, not a modification of a standard NASTRAN solution. Like 
nlcfs, both quadratic a jk and cubic b jkl stiffness coefficients are computed in nlcfs.bdf, but only 

the cubic terms are stored in nlcfS.dat because a zero mean response is assumed. The nonlinear 
stiffness coefficients are also stored in the ASCII file nlcfS.dat, but in the INPUTT4 DMAP 
matrix input format. 

As in the FORTRAN implementation, this portion of the analysis only needs to be performed once for a 
given structural configuration. Similarly, proper functioning of the entire analysis procedure should be 
checked using the method indicated in Section 3.2. 

The UNIX script file doitDl.bat (see Appendix A) runs the first part of the analysis as outlined above in 
an automated fashion. 
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Figure 15: Flow-chart for calculation of nonlinear stiffness coefficients in the DMAP implementation. 
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4.3. Equivalent Linearization Procedure 

Like the FORTRAN implementation, both force and strain energy error minimization approaches follow 
the same sequence of analyses, as shown in Figure 16. Each operation is separated by a dashed line and is 
run in the sequence shown from top to bottom. The various operations are described below. 


OPTIONAL 



OPTIONAL 



Figure 16: Flow chart for calculation of total equivalent stiffness matrix and RMS displacements in the 

DMAP implementation. 

• gensffD - The first operation is optional and is used to develop the spectral density matrix of the 
loading in modal coordinates. It is analogous to the program gensff in the FORTRAN 
implementation, except that in the DMAP implementation, the spectral density matrices are 
independent of frequency, i.e. only flat spectra may be specified. The program gensffD may be 
used to compute the modal spectral density matrix for a couple of special cases - a uniformly 
distributed acoustic pressure (which may also be used to specify an inertial load) and point forces 
and moments. In each case, the loadings are fully correlated giving a real-valued spectral density 
matrix in physical coordinates. The mass-normalized selected set of g-size eigenvectors are read 
from egveCa.dat and information about the loading is read from fixeD.sff. The computed modal 
spectral density matrix is written to fixeD.den (as described in Section 4.1.2). 

o The parameter file fixeD.sff specifies the following information in free format: 

LINE 1 : Load option (LO ADOPT) 

1 = uniformly distributed pressure with uniform mesh 

2 = point force(s) and moment(s) 

If LOADOPT = 1 

LINE 2: PSD, AREA, NOD, NODB 
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PSD - Spectrum level of uniformly distributed pressure in units 2 /Hz. This level 
is a single-sided value and is the same as would be used on the TABRND1 card 
for sol 1 1 1 analysis. Note that gensffD converts the load internally to a double- 
sided value in units 2 /rad/sec. 

AREA - Surface area of structure 

NOD - Total number of nodes in the model 

NODB - Number of constrained (boundary) nodes 

Forces at the grid points are obtained by multiplying the pressure level by the 
area and distributing these over the mesh, hence the need for a uniform mesh. If 
a non-uniform mesh is used, a pressure loading may be specified using 
LOADOPT = 2 with forces input at each grid point. 

Note that in this case of a uniformly distributed pressure on a uniformly 
distributed mesh, the forces for each grid point are identical and act in the z- 
direction. Hence, all non-zero auto- and cross-spectral density values in S (f are 
equal. Further, each column is identical. These properties speed the 
computation of S ^ . 

If LOADOPT = 2 
LINE 2: NLOAD 

NLOAD - Number of fully correlated point forces and moments. 

LINE 3 (repeated NLOAD times): NODF, PSDF, DOF1-DOF6 
NODF - Grid point at which point force is applied 

PSDF - Spectrum level of force in units 2 /Hz. This level is a single-sided value 
and is the same as would be used on the TABRND1 card for sol 1 1 1 analysis. 
Note that gensffD converts the load internally to a double-sided value in 
units 2 /rad/sec. 

DOF1-DOF6 - Multiplier for each of six degree-of-freedoms at each grid point. 
An arbitrarily oriented force or moment may be applied in this fashion. 

Note that in this case, the forces for each grid point are not identical. Therefore, 
the non-zero auto- and cross-spectral density values in Sff are unequal and each 

column is not identical. These properties make the computation of S ff more 
intensive. 

Note that in general, the cross-spectral densities of the input loading are required to 
specify S (f . However, since both load options require fully correlated (in phase) 

loadings, the cross-spectral density (CSD) terms may be obtained from the auto-spectral 
density (PSD) from the following relationship: 

CKD,=±^(«D,)(PSD,) 

This allows for the specification of only the PSD values, and not both the PSD and CSD 
values, as indicated above. Note that gensffD as currently implemented only takes the 
positive square root. Therefore, the specified loadings must have the same direction. 
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Failure to do so will result in an inconsistency in the sign of the diagonal and off-diagonal 
terms in S ff . 

The following is an example of a uniformly distributed pressure of 4096 PaVHz: 

1 

4096 0.0903224 609 96 

The following is an example of a point force of 5 N7Hz in the z-direction at grid point 
305 and a point force of 2 N 2 /Hz in the x-direction at grid point 220: 

2 


2 

305 5.0 

220 1 - 0 


0 . 0 
2 . 0 


0 . 0 
0 . 0 


1.0 
0 . 0 


0 . 0 
0 . 0 


0 . 0 
0 . 0 


0 . 0 
0 . 0 


• prepelD - This FORTRAN program is the analogue to prepel in the FORTRAN implementation. 

It reads the NASTRAN f06 file fixeD103_eg.f06 produced in the second step of the previous 
procedure (see Section 4.2) and writes the selected natural frequencies to the file fixeD.stt. Here, 
fixeD.stf is an ASCII file with the INPUTT4 DMAP matrix input format. 

. The third operation is the DMAP analogue of the FORTRAN programs eqlf and eqb. This 
NASTRAN run performs the equivalent linearization procedure through the standalone DMA! 
alters eqlj.bdf and eql_s.bdf, for the force error and strain energy error minimization 
approaches, respectively. Each of these alters produce the following output: 

o st f.dat and st_s.dat contain the total equivalent stiffness matrix in modal coordinates for 
the force and strain energy error minimization approaches, respectively. These are ASCII 
files with the INPUTT4 DMAP matrix input format. 

o rms f.dat and rms_s.dat contain the root-mean-square displacements in physical 

coordinates for all degrees of freedom. These are ASCII files with the INPUTT4 DMAP 
matrix input format. 

o dcov_f.dat and dcov_s.dat contain the cross-covariance matrix of modal displacements. 
These are ASCII files with the INPUTT4 DMAP matrix input format. 

• The fourth (optional) operation has no FORTRAN analogue. It should be performed only if 

graphical display of results in PATRAN is desired. FORTRAN programs postelfD and postelsD 
read the files rms f.dat and rms_s.dat, respectively, and separate the translational and rotationa 
degrees of freedom. The RMS translational displacements for all grid points are written to either 
rms f.nod or rms_s.nod, depending on the error minimization procedure used. The RM. 
rotational displacements for all grid points are written to either rms JR.nod or rms_sR.nod 
depending on the error minimization procedure used. These ASCII results files may be importe 
from within PATRAN for plotting, etc. 

There are a couple of comments worth noting about the analysis described in this section. The first has to 
do with creation of the files eql_f.f04 and eql_s.f04 when running eqlj.bdf and eq1_s.bdf respectively. 
These files can be very large due to the amount of information written by NASTRAN. Unfortunately, a 
means of suppressing creation of these files or minimizing the information written to them has not been 
identified. Although the files can be deleted immediately following successful execution of the run 
ample disk space must be available to successfully complete the run. This problem is most acute w en 
running the strain energy error minimization analysis, as the number of iterations for convergence at 
highly nonlinear levels can be in the hundreds or thousands (compared with ten or less typical of the force 
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error minimization approach). Consequently, large temporary disk space is recommended when running 
the strain energy error minimization analysis. 

The second item worth mentioning relates to the calculation of RMS displacements in physical 
coordinates when running eql_f.bdf and eql_s.bdf. Following calculation of the total equivalent stiffness 
matrix and cross-covariance matrix in modal coordinates, the RMS displacements in physical coordinates 
are recovered and stored in the files rms_f.dat and rms_s.dat. This process can take a significant amount 
of time as presently implemented due to the amount of disk I/O that occurs, especially for large models. 

If the intention is to perform a post-processing analysis to determine displacement PSDs and/or 
stress/strain information, as outlined in Section 5, a more efficient means of getting RMS displacements 
in physical coordinates is to request this information as part of that analysis. In this case, the value of 
RMSDISP should be changed from 1 to 0 in eql _f.bdf and eql_s.bdf to suppress computation of this data. 

The UNIX script files doitD2_f.bat and doitD2_s.bat (see Appendix A) run the second part of the 
analysis as outlined above in an automated fashion. 


4.3.1. Convergence 

The reader should refer to Section 3.3. 1 for a discussion on convergence of the equivalent linearization 
procedure. 
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5. Post-Processing for Displacement, Stress and Strain PSDs 

An optional post-processing procedure may be performed to compute RMS displacements, displacement 
PSDs , RMS stress/strain, and stress/strain PSDs. It applies to both FORTRAN and DMAP 
implementations. A flowchart of the implementation is shown in Figure 17. Each operation is separated 
by a dashed line and is run in the sequence shown from top to bottom. The various operations are 
described below. 



Figure 17: Flow chart for post-processing. 

• The first operation simply copies the total equivalent stiffness matrix from either the force or 
strain energy error minimization approaches to the file st.dat. This allows stiffness matrices from 
both approaches to coexist in the same directory. 

• The second operation performs a post-processing analysis using the total equivalent stiffness 
matrix, obtained from the equivalent linearization analysis, in place of the linear stiffness matrix 
in a NASTRAN random analysis (solution 111). 

Note that since the total equivalent stiffness matrix in modal coordinates is symmetric, but not (in 
general) diagonal, the modal system of equations to be solved is coupled. In order to reduce the 
computational effort, the option is provided to analyze an uncoupled system, which is formed 
through an additional transformation process. The uncoupled system has a different total 
equivalent stiffness matrix (now diagonal and in transformed coordinates), but utilizes the same 
diagonal damping matrix and unit mass matrix as the coupled system. Analysis of the uncoupled 
system is significantly faster than the coupled system, but is not as accurate. It is recommended 
that the root-mean-square displacements from the post-processing analysis be compared with 
those obtained from the equivalent linearization analysis. In all circumstances, post-processing 
results from the coupled system should be comparable to values obtained from the equivalent 
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linearization analysis. Only small differences, on the order of a few percent, should be found (see 
Section 6.4 for example), since the loading is specified equivalently but differently between the 
equivalent linearization and post-processing analyses. Post-processing results from the 
uncoupled system may or may not be comparable to values obtained from the equivalent 
linearization analysis. Thus, if results from the uncoupled analysis significantly differ from 
results from the equivalent linearization analysis, it is recommended that the coupled system be 
used for post-processing. The choice of coupled versus uncoupled system is made in the file 
fixeDl llrd.bdf (see Appendix D), through selection of the appropriate DMAP alter as shown in 
Figure 17. 

If root-mean-square strain results are to be recovered, the alters strain700.v or strain705.v 
should be utilized. For NASTRAN version 70.0, the alter strain700.v should be used. For 
NASTRAN versions 70.5 and 70.7, the alter strain705.v should be used. If root-mean-square 
stress results are to be recovered, neither alter should be utilized. The choice of which alter to use 
(if any) is made in the file fixeDl 11 rd.bdf (see Appendix D). The file Strain-Recovery- 
Notes.pdf (see Appendix A) documents other commands required in the case control section for 
the different NASTRAN versions. Note that stress and strain results are recovered using linear 
strain/displacement relationships. This approach appears to be acceptable for the beam and plate 
analyses conducted to date. 

Because the post-processing analysis is performed on the equivalent linear system, all PSDs 
appear linear in character. In other words, at nonlinear response levels, the PSDs exhibit neither 
peak broadening nor the same shifts in resonant frequencies as the truly nonlinear system. 

Shifted resonant frequencies of the equivalent linear system are not directly reported since a 
modified solution 103, using the total equivalent stiffness matrix in place of the linear stiffness 
matrix, is not presently available. The shifted resonant frequencies may be determined by 
scanning the punch file fixeDl 1 l_rd.pch for peaks in the displacement PSD (when flat spectra 
are applied). Note that it may be necessary to run the post-processing analysis more than once so 
that the frequency ranges specified on the TABDMP1 card are centered about the shifted peaks of 
the equivalent linear system. The frequency range for each critical damping value is typically set 
to bracket the corresponding resonant peak at mid-span. As indicated in Section 3. 1 .3, the 
relationship between percent critical damping and modal damping values should remain based on 
the linear eigenvalues. 

The UNIX script files doitD3_f.bat and doitD3_s.bat (see Appendix A) run the post-processing part 

of the analysis as outlined above in an automated fashion. 
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6. Example Problem 


A simply support aluminum plate under normal acoustic excitation was chosen to illustrate the execution 
of ELSTEP and to generate sample results which may be used subsequently to verify proper execution 
with future NASTRAN and PATRAN releases. Consequently, many results are reported to a large 

number of significant digits. 

The rectangular plate chosen measured 0.254m x 0.3556m x 0.00102m. All srdes were simply supported. 
Material properties were: 


E = 7.3x 10 10 Pa 
V = 0.3 


giving a flexural rigidity of 


p = 2763 kg/m 3 


D = 


Elf 

12(1-F 2 ) 


= 7.094 


kg-m 3 

i 

s" 


Natural frequencies may be computed from the following relation [ 1 5] 

[5 


mK 

a ) 


Mt 


giving 


an, = 58.34 Hz 
co n = 216.01 Hz 


for the first two symmetric modes. 

A finite element model of this structure was built in PATRAN having 560 CQUAD4 elements measuring 
0.0127m square. The first two symmetric modes (modes 1 and 4) of this model had natural 
58 38 and 217 27 Hz These two modes were selected as participating modes in the ELS y ^ < - 

per the example files in Sections 3 . 1 . 1 and 4. 1 . 1 . Modal damping was specified as per the example in 
Sections 3. 1 .3 and 4. 1 .3, giving percent critical damping values of 2% and 0.54 /o for modes and 4, 

respectively. 

A spatially uniform distributed normal acoustic pressure with RMS value of 2048 Pa, and a single^ided 
spectrum level of 4096 Pa’/Hz over a 0- 1 024 Hz bandwidth was applied to the plate. Using the example 
for a uniformly distnbuied pressure (described in Sections 3.3 and 4 or the I programs gensff 
and gensffD, the modal spectral density matrix generated is as shown in Sections 3.1 .2 a ... 
that this load level gives rise to a highly nonlinear response. 


6.1. Nonlinear Stiffness Results 

The cubic nonlinear stiffness coefficients were computed and written to nlcfS.dat. For the FORTRAN 
implementation, nlcfS.dat is written as. 


4 . 1355686358055371E+12-1 . 6827938202459033E+12 
2 .3630088637711590E+13 
4 . 0490613582441501E+07 
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-5 . 60514 5 07 03 7 9141 8E+11 2 . 3 62 8 886 15 86854 77E+ 13 
3 . 58870215279804 59E+ 08 
7. 047404098343 54 22 E+ 13 

For the DMAP implementation, nlcfS.dat is written as: 

441 2 CUB IP, 3E23 . 16 

112 

4 . 1355688017920000E+12-5 . 6051453132800000 E+ 11 

3 14 

-1 . 6827939553280000E+12 2 . 3628886310912000E+13 2 . 3630090076160000E+13 
3 . 5887843200000000E+08 

4 3 2 

4 . 04 9 0628 00 0000 000E+ 07 7 . 04 74 046373 8 8 800 0E+ 13 
In less cryptic form, the stiffness terms for the first selected mode may be written as: 

= +4.1356x1 O' 2 
b' m = - 1.6828 xlO 12 
b\ | 22 = +2.3630 xlO 13 

^222 =+4 0491Xl()7 

and for the second selected mode may be written as: 

bf u =-5.6051x10" 
bf ]2 = + 2.3629 xlO 13 
bf 22 = +3.5888x10* 

^222 =+ 7 0474x1 °' 3 

It is seen that the FORTRAN and DMAP implementations give essentially the same nonlinear stiffness 
coefficients. 

6.2. Equivalent Linearization Results - Force Error Minimization 

Results from the equivalent linearization analysis using the force error minimization approach were 
generated using the example parameter files provided in Sections 3. 1 .4 and 4. 1 .4. For both FORTRAN 
and DMAP implementations, convergence was achieved after ten iterations. 

From the FORTRAN implementation, the covariance matrix of modal displacements is given as: 

'3 .3063497336422452 x 1 0‘ 7 7.792070656260 1 622 x 1 0‘ 9 " 

7.792070656260 1 688 x 1 0‘ 9 1 . 1 782 1 92964404840 x 1 0' 8 
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For the DMAP implementation, the nearly identical covariance matrix is given as: 

3.306347272 1 539035 x 10 7 7.7920762241445690 x 10 9 

7.792076224 1 445690 x 1 0 9 1.1 782147220664467 x 1 0 8 

The total equivalent stiffness matrix in modal coordinates from the FORTRAN implementation is given 
as: 

" 4.5034794166661 100x1 0 6 -1. 843 5794253364013 x 10 5 

-1.843579425336401 3x 10 5 1.20473868231 1 3656x1 0 7 

The DMAP implementation gives very similar values of 

’ 4.49761 86005045557 xlO 6 -1 .8586848433203337 x 10 5 
- 1 .8586848 173153755x1 0 5 1 .20953 1 057076 1 03 1 x 1 0 7 

The RMS displacement at the plate center (node 305) from both FORTRAN and DMAP implementations 
is 2.27352 x 1 0' 3 m. A linear analysis at this level would produce an RMS center displacement of 
1.31952X 10~ 2 m or w miX /h of 12.9. A nonlinear analysis is typically deemed necessary for values of 

w max //? in excess of 0.5 for plate analysis. 

6.3. Equivalent Linearization Results - Strain Energy Error Minimization 

Results from the equivalent linearization analysis using the strain energy error minimization approach 
were generated and are provided below. Due to the high degree of nonlinearity, values of a and (5 of 
0.0005 and 0.9995 were used to achieve a monotonically decreasing convergence error. Convergence 
was achieved after 4938 iterations. Because of the large number of iterations, only the more 
computationally efficient FORTRAN implementation was invoked. 

The covariance matrix of modal displacements is given as: 

3.4823285014925753 xlO' 7 7.3445452335136257 x 1 O' 9 
7.3445452335 1 30459 x 1 O' 9 2.590 1 93 1 870425043 x 1 O' 8 

and the total equivalent stiffness matrix in modal coordinates is given as: 

' 4.2663 1 54 1 833487 1 5 x 1 0 6 -2.748647453 1 34 1 883 x 1 0 4 " 

-2.748647453 1 34 1 883 x 1 0 4 5.5 145054668358499 x 1 0 6 

The RMS displacement at the plate center is 2.38505 x 10' 3 m. The strain energy error minimization 
approach is seen to give more conservative results than the force energy minimization results (typically 
5% higher), albeit at a somewhat higher computational expense due to the large number of iterations 
required for convergence. The strain energy minimization results have been found to compare more 
favorably with fully nonlinear analyses than force error minimization results [10], 

6.4. Post-Processing Results 

The total equivalent stiffness matrices from the force and strain energy minimization equivalent 
linearization approaches were used to post-process RMS displacements, displacement PSDs, and stress 
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PSDs. The following uncoupled post-processing analysis results are derived from the total equivalent 
stiffness matrices obtained from the FORTRAN implementation. 


As previously suggested, RMS displacements were computed to serve as a consistency check with the 
equivalent linearization results. These are summarized in Table 1 and are seen to fall within the typical 
range of differences. 


Center Displacement (m) Center Displacement (m) 
EL Approach f rom el from Post-Processing 


% Difference 
in Center 
Displacement 


Force Error 2.2735 x 10' 3 2.2463 x 10' 3 

Minimization 


1.2 


Strain Energy Error 2.3851 xlO 3 2.3200 x 1 O' 3 2.7 

Minimization 

Table 1 : Comparison of RMS displacements from equivalent linearization and post-processing. 


A plot of the center displacement PSDs for the two equivalent linearization approaches is shown in Figure 
18. These results show that the resonant frequencies of the equivalent linear system using force error 
minimization shift higher to a greater degree than do those associated with the strain energy error 
minimization approach. This is as expected since the force error minimization approach results in a 
stiffer equivalent linear system. Recall that since the two equivalent linearization analyses only minimize 
the error between the linear and equivalent linear systems, only the RMS displacements should be 
comparable, and not the shape of their PSDs. 

A plot of the maximum element stress PSD (taken for plate center element 266 across the short 
dimension) is shown in Figure 19. Similar trends as those above are also noted in the stress PSDs. The 
RMS levels are comparable - 15.92 MPa for the force error minimization versus 16.73 MPa for the strain 
energy minimization. For comparison, a linear analysis predicts 95.56 MPa. 
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Frequency (Hz) 


Figure 18: Comparison of displacement PSDs using linear and equivalent linear analyses. 



Frequency (Hz) 

Figure 19: Comparison of stress PSDs using linear and equivalent linear analyses. 
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7. Summary 


An equivalent linearization procedure using nonlinear stiffness evaluation for analysis of geometrically 

nonlinear random vibrations has been implemented in two forms. Both ' ^ r ^ r ANTmn le menfat . on 
commercial finite element program MSC/NASTRAN as their core. The FORTRAN implementation 

calculates the nonlinear stiffness terms and performs the equivalent 

NASTRAN. The DMAP implementation performs these opeiations within • P 

normal modes (solution 103), linear static (solution 101) and nonlinear static (solution 106) analyses 
thtTcalculatlori of the nonlinear s>,ffn=s s lerms. Optional post-procasstng is performed by subs trtutrng the 
equivalent stiffness matrix in place of the linear stiffness matrix in a random (solution 1 1 1) analysis. The 
following items are noteworthy: 

• Because of the way NASTRAN solvers are implemented, a second model database is required 
which adds a dummy element and force. These are required to solve the linear and nonlinear 
inverse static problems for calculation of the nonlinear stiffness terms. 

• Both FORTRAN and DMAP implementations allow specification of real-valued spectral density 
matrices of the applied loading in physical coordinates. This permits only correlated inputs sue 
as normal acoustic pressure, inertial loadings, and one or more correlated point forces or 
moments. Loadings giving rise to complex-valued spectral density terms, such as grazing 
incidence acoustic loading, are not permitted. 


• The FORTRAN implementation allows for a frequency dependent spectral density matrix of the 
applied loading. The DMAP version allows for only a flat input spectrum. 

• Both implementations allow force and strain-energy error minimization approaches in the 
equivalent linearization procedure. 

. Results are comparable between the FORTRAN and DMAP implementations. 

• Convergence is typically on the order of ten iterations for the force error minimization approach 
and can be in the thousands of iterations for the strain-energy error minimization approach. 

• Within the equivalent linearization analysis, the nonlinear solution space should be explored 
using different weightings to determine whether multiple solutions exist. 

• Results from the strain-energy error minimization approach have been previously shown to be 
closer to a fully nonlinear analysis than are those from the force error minimization approach. 

. Results from the strain-energy error minimization approach are more conservative (typically 5% 
higher) than are those from the force error minimization approach. 

• Post-processing results based upon the coupled equivalent linear system are more accurate than 
those based upon an uncoupled system, but are obtained with greater computational effort. 


• The FORTRAN implementation is more computationally efficient than the DMAP 
implementation. 
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Appendix A Software Distribution and Installation 


Distribution 

Requests for the ELSTEP software distribution should be submitted in writing to the following address: 

Dr. Stephen A. Rizzi 

NASA Langley Research Center 

Mail Stop 463 

Hampton, VA 2368 1-2199 

Email: s.a.rizzi@larc.nas a.gov 

The request should include the requester’s complete name, address, telephone and fax numbers, email 
address, and citizenship. The software will be distributed in an appropriate form following completion 
and acceptance of a software usage agreement. 


Installation 

The ELSTEP source code is organized in the following tree structure: 

README.TXT 

dmap 

README.TXT 

dmlOlD.v 

dm_103D.v 

dm_106D.v 

doitDl.bat 

doitD2_f.bat 

doitD2_s.bat 

doitD_f.bat 

doitD_s.bat 

eqljlbdf 

eql_s.bdf 

gensffD.f 

makeall.bat 

makefile 

nlcfs.bdf 

param.h 

postelfD.f 

postelsD.f 

preplOlD.f 
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prepl03D.f 

prepl06D.f 

prepelD.f 

doc 

README.TXT 

NASA-TM-2001-210838.pdf (this document) 
MSC_Conf_6-99.pdf 
NASA-2000-7icrasd-sar.pdf 
NASA-aiaa-99-1 376.pdf 
Strain-Recovery-Notes.pdf 

example 

README.TXT 

dmap 

fixeD.dam 
fixeD.den 
j fixeD.mod 
fixeD.par 
fixeD.sff 

fortran 

README.TXT 
fixed.dam 
fixed.den 
fixed.mod 
fixed. parf 
fixed. pars 
fixed. sff 

patran 

ssplatel_eg.bdf 

ssplatel_eg.db.jou 

ssplatel_rd.bdf 

ssplatel_st.bdf 

ssplatel_st.db.jou 

post 

README.TXT 

fixeD.mod 
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fixeDl 1 Ird.bdfs 
fixeDlll rd.bdff 


README.TXT 

cdinvt.f 

drrMOl.v 

dm_103.v 

dm_106.v 

doitl.bat 

doit2__f.bat 

doit2_s.bat 

eqlf.f 

eqls.f 

e qiii v _st,f 

gendisp.f 

gensff.f 

inv.f 

makeall.bat 
makefile 
nlcfs.f 
param.h 
preplOl.f 
prep 1 03. f 
prep 1 06. f 
prepel.f 


README.TXT 
dm_l 1 lu.v 
dm_l 1 lc.v 
strain700.v 
strain705.v 
doitD3_f.bat 
doitD3 s.bat 
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The dmap directory off the main tree contains all the files required to run the DMAP inylenMmtadon of 
the analysis All FORTRAN utility programs contained in this directory are platform independen 
should be compiled with a FORTRAN 90 compiler. The batch file makefile.bat will compile all 
FORTRAN programs using the UNIX “make” utility. Prior to compilation, the values fortbem^imum 
number of elements in the model (MAXEL), the maximum number of degrees of freedomm themodel 
(NSX) the maximum number of eigenvectors specified in the normal modes analysis (NEIGX), and the 
maxfmum number of selected eigenvectors (NSELX) should be set sufficiently large in the file param.h. 
The DMAP implementation uses string-based Alters of NASTRAN version 70.0.0. It is therefore 
expected to be upward compatible with future versions of NASTRAN. Following compilation, the 
should be copied to the directory in which the analysis is to be performed. 

The fortran directory off the main tree contains all the files required to run the FORTRAN 
implementation of the analysis. All FORTRAN programs contained in this directory are platform 
independent and should be compiled with a FORTRAN 90 compiler. The batch fi e makefile.bat will 
compile all FORTRAN programs using the UNIX “make” utility. Prior to compilation, the values for the 
maximum number of elements in the model (MAXEL), the maximum number of degrees of freedom in 
the model (NSX), the maximum number of eigenvectors specified in the normal modes analysis 
(NEIGX) and the maximum number of selected eigenvectors (NSELX) should be set sufficiently large in 
ihe file param.h. Following compilation, the files should be copied to the directory in which the analysis 

is to be performed. 

Note that since some of the files created in the DMAP and FORTRAN implementations have the same 
name, but differing content, the two implementations should be run in separate directories. 

The post directory off the main tree contains all the files required to run the post processing analysis for 
either DMAP or FORTRAN implementations. These files should be copied to the directory in which the 
DMAP or FORTRAN analyses are performed. No compilation is required. 

The doc directory off the main tree contains an Adobe Acrobat version of this document as well as 
Acrobat versions of references [10-12], 

The example directory off the main tree contains the files required to run the simply support plate 
example problem used in this document. The dmap, fortran, and post sub-directones contain The mpirt 
files required to run the DMAP, FORTRAN and post-processing analyses, respectively. The files in these 
directories should be copied to the appropriate analysis directory. The patran sub-directory contains 
PATRAN bulk data files and database journal files for the simply supported plate model men miming 
the example problem, two bulk data files (ssplatel_eg.bdf and ssplatel_st.bdf) should be copied to the 
analysis directory. The third bulk data file, ssplatel_rd.bdf, has already modified for post-processing as 
indicated in Section 2.3. The modified versions of this file are m the post sub-directory as 
fixeDl 1 1 rd.bdfs and fixeDl 1 l_rd.bdff for the strain-energy and force error nunimrzation approaches. 
The only difference between the two files are the TABDMP 1 card which indicates slightly different 
damping bounds to reflect the different shifts in fundamental frequency associated with the two Afferent 
error minimization approaches. The appropriate one should be copied to the file fixeDl 11 jrd.bdf in 
order to run the post-processing analysis. The PATRAN databases themselves may be re-created (if 
desired) using either the bulk data files or the database journal files. 
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Appendix B Listing of Sample Normal Modes Analysis Bulk Data Fde 


$ NASTRAN input file created by the MSC MSC/NASTRAN input file 
$ translator ( MSC/PATRAN Version 8.0 ) on October 27, 20 

$18:06:02. 

ASSIGN OUTPUT2 = ' ssplatel_eg . op2 ' , UNIT = 12 

$ Direct Text Input for File Management Section 
$ Normal Modes Analysis, Database 
SOL 103 

TIME 600 . . 

$ Direct Text Input for Executive Control 

CEND 

SEALL - ALL 
SUPER — ALL 

TITLE = MSC/NASTRAN job created on 27 -Oct- 00 at 15:49:00 
ECHO = NONE 
MAXLINES - 999999999 

$ Direct Text Input for Global Case Control Data 
SUBCASE 1 

$ Subcase name : Eigen 
SUBTITLE= Eigen 
METHOD = 1 
SPC = 2 

VECTOR (SORT 1 , REAL) =ALL 
$ Direct Text Input for this Subcase 
BEGIN BULK 
PARAM POST - 1 

PARAM PATVER 3 . 

PARAM AUTOS PC YES 

PARAM COUPMASS 1 
PARAM K6ROT 0 . 

PARAM WTMASS 1. 

PARAM , NOCOMPS , - 1 
PARAM PRTMAXIM YES 

EIGRL 1 10 0 

$ Direct Text Input for Bulk Data 
$ Elements and Element Properties for region 
1 .00102 1 

will be imported as 
1 1 2 


MASS 


PSHELL 

$ Pset : 

CQUAD4 

CQUAD4 

CQUAD4 

CQUAD4 


1 

n el prop" 
1 
2 

3 

4 


"pshell . 1" 

31 

32 

33 

34 


elprop 

1 

30 

31 

32 

33 


576 

577 

578 

579 


C QUAD 4 557 1 

CQUAD4 558 1 

CQUAD4 559 1 

CQUAD4 560 1 

$ Referenced Material Records 
$ Material Record : alum 
$ Description of Material : Date 
MAT1* 1 7.3+10 

* A 2763 . 

$ Nodes of the Entir 
GRID 1 

GRID 2 

GRID 3 

GRID 4 


577 

578 

579 

580 


606 

607 

608 
609 


27 -Oct -00 


0 . 

0 . 

0 

. 0127 

0 . 

0 

. 0254 

0 . 

0 

. 0381 

0 . 

0 


605 

606 

607 

608 


Time : 15:41:56 

. 3 


GRID 

GRID 

GRID 

GRID 


606 

607 

608 
609 


$ Loads for Load Case 


. 3175 
. 3302 
. 3429 
.3556 
Eigen 


. 254 
. 254 
.254 
.254 
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SPCADD 2 

$ Displacement 
SPC1 1 

+ B 175 

+ C 407 

$ Displacement 
SPC1 3 

+ D 203 

+ E 435 

$ Displacement 
SPC1 4 

$ Displacement 
SPC1 5 


1 3 

Constraints of 
134 1 

204 233 

436 465 

Constraints of 
134 29 

232 261 

464 493 

Constraints of 
235 1 

Constraints of 
235 581 


4 5 


Load Set : 

: xOsimple 




30 

59 88 

117 

146 

+ 

262 

291 320 

349 

378 

+ 

494 

523 552 

581 



Load Set 

: xl simple 




58 

87 116 

145 

174 

+ 

290 

319 348 

377 

406 

+ 

522 

551 580 

609 



Load Set 

: yOsimple 




THRU 

29 




Load Set 

: ylsimple 




THRU 

609 





$ Referenced Coordinate Frames 
ENDDATA 15845flb 


B 

C 


D 

E 
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Appendix C Listing of Sample Nonlinear Static Analysis Bulk Data File 


$ NASTRAN input file created by the MSC MSC/NASTRAN input file 
$ translator { MSC/PATRAN Version 8 . 0 ) on October 27, 2000 at 
$ 16:51:35. 

ASSIGN OUT PUT 2 = ' ssplatel_st . Op2 ' , UNIT = 12 
$ Direct Text Input for File Management Section 
$ Nonlinear Static Analysis, Database 
SOL 106 
TIME 600 

$ Direct Text Input for Executive Control 
CEND 

SEALL = ALL 
SUPER - ALL 

TITLE = MSC/NASTRAN job created on 27-Oct-OO at 16:47:55 
ECHO = NONE 
MAXLINES = 999999999 

$ Direct Text Input for Global Case Control Data 
SUBCASE 1 

$ Subcase name Default 
SUBTITLE=Def aul t 


NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1 , REAL) =ALL 
$ Direct Text Input for this Subcase 


BEGIN 

BULK 


PARAM 

POST 

- 1 

PARAM 

PATVER 

3 . 

PARAM 

AUTOS PC 

NO 

PARAM 

COUPMASS 

1 

PARAM 

K6ROT 

100 

PARAM 

WTMASS 

1 . 

PARAM 

LGDISP 

2 

PARAM , NOCOMPS , - 1 


PARAM 

PRTMAXIM 

YES 


NLPARM 1 1 AUTO 

* A 1.+13 

+ B 


1 50 

1 . +13 


PW 


$ Direct Text Input for Bulk Data 
$ Elements and Element Properties 


for region 


PSHELL 1 1 

$ Pset: "elprop" will be 

C QUAD 4 
CQUAD4 
CQUAD4 
CQUAD4 


.00102 1 
imported as: 


elprop 

1 


"pshel 1 . 1 " 


1 

1 

1 

2 

31 

30 

2 

1 

2 

3 

32 

31 

3 

1 

3 

4 

33 

32 

4 

1 

4 

5 

34 

33 

557 

1 

576 

577 

606 

605 

558 

1 

577 

578 

607 

606 

559 

1 

578 

579 

608 

607 

560 

1 

579 

580 

609 

608 


Properties for 
.01 

be imported as : 
1 610 


CQUAD4 
CQUAD4 
CQUAD4 
CQUAD4 

$ Elements and Element 
PROD 2 1 

$ Pset: "rodprop" will 

CROD 561 2 1 

$ Referenced Material Records 
$ Material Record : alum 
$ Description of Material : Date 
MAT1 * 1 7.3+10 

* C 2763. 

$ Nodes of the Entire Model 
GRID 1 0* 

GRID 2 .0127 

GRID 3 .0254 

GRID 4 .0381 


region : rodprop 
"prod . 2 " 


27 -Oct - 00 


Time : 
. 3 


NO 


16 : 28 : 40 


A 

B 


C 


c-i 



GRID 

606 


.3175 

. 254 

GRID 

607 


.3302 

. 254 

GRID 

608 


. 3429 

.254 

GRID 

609 


.3556 

.254 

GRID 

610 


- . 05 

0 . 

$ Loads 

for Load Case : 

Default 

SPCADD 

2 

1 

3 

4 

LOAD 

2 

1 . 

1 . 

1 

$ Displ 

acement 

Constraints of 

Load Set 

SPC1 

1 

134 

1 

30 

+ 

D 175 

204 

233 

262 

+ 

E 407 

436 

465 

494 

$ Displacement 

Constraints of 

Load Set 

SPC1 

3 

134 

29 

58 

+ 

F 203 

232 

261 

290 

+ 

G 435 

464 

493 

522 

$ Displacement 

Constraints of 

Load Set 

SPCI 

4 

235 

1 

THRU 

$ Displacement 

Constraints of 

Load Set 

SPCI 

5 

235 

581 

THRU 

$ Displacement 

Constraints of 

Load Set 

SPCI 

6 

123456 

1 

THRU 

$ Nodal 

Forces 

of Load 

Set : dummy_f 

FORCE 

1 

610 

0 

12 . 


$ Referenced Coordinate Frames 
ENDDATA dbd9fe3d 


0 . 

0 . 

0 . 

0 . 

0 . 

5 6 


xOsimple 
59 88 

117 

146 

291 

320 

349 

378 

523 

552 

581 


xls impl e 
87 116 

145 

174 

319 

348 

377 

406 

551 

580 

609 


yOsimple 

29 

ylsimpl e 
609 

const rained 
609 




1 . 0 . 0 . 


F 

G 
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Appendix D Listing of Sample Random Analysis Bulk Data File 


ASSIGN OUTPUT 2 

ASSIGN INPUTT4 
ASSIGN INPUTT4 
ASSIGN INPUTT4 


$ NASTRAN input file created by the MSC MSC/NASTRAN input file 
$ translator ( MSC/PATRAN Version 8.0 ) on October 
e 17-57:54. 

rfixeDlll rd . op2 ' , UNIT =12 

1 f ixeD .mod' , UNIT =65, FORM = FORMATTED 
•st dat', UNIT = 82, FORM = FORMATTED 
•egveH.dat', UNIT = 83, FORM = UNFORMATTED 

$ Direct Text Input for File Management Section 
$ Frequency Response Analysis, Modal Formu ation, 

SOL 111 
TIME 600 

$ 

$ For an uncoupled system, 

INCLUDE ' dmlllu.v' 

$ OR, for a coupled system 
INCLUDE ' dm_lllc . v ’ 


Database 


include the following 
, include the following 


For strain recovery, include one of the following: 
For NASTRAN version 70.0 - 
INCLUDE ' strain700 . v' 

For NASTRAN versions 70.5 or 70.7 - 
INCLUDE ' strain705 . v* 


$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ Direct Text Input for Executive Control 
CEND 

$ Elements for group : centerelem 
SET 1 = 266,267,294,295 
$ Nodes for group : center_nodes 

SET 2 = 275,276,277,304,305,306,333,334,335 

SEALL = ALL 

TITLE = MSC/NASTRAN job created on 27-Oct-OO at 15:49:00 

ECHO = NONE 
MAXLINES = 999999999 

$ Direct Text Input for Global Case Control Data 
SUBCASE 1 

$ Subcase name : random 
SUBTITLE= random 


METHOD = 1 
FREQUENCY = 2 
SPC = 2 
LOAD SET = 1 
DLOAD = 2 

DISPLACEMENT {SORT1 , REAL) =2 
STRESS ( SORT1 , REAL, VONM I SES , CENTER) =1 
$ STRAIN (SORT1, PLOT, REAL , FIBER) =1 
$ Direct Text Input for this Subcase 

SDAMP ING= 4 0 0 
RANDOM= 1 
OUTPUT (XYPLOT ) 

PLOTTER NAST 

$ Following lines specific to desired output 
XYPUNCH DISP PSDF / 305 (T3) 


XYPUNCH 

STRESS PSDF / 

266 (3) 

XYPUNCH 

STRESS PSDF / 

266 (5) 

XYPUNCH 

STRESS PSDF / 

266 (7) 

BEGIN BULK 



PARAM 

POST 

-1 


PARAM 

PATVER 

3 . 


PARAM 

AUTOS PC 

YES 


PARAM 

COUPMASS 

1 


PARAM 

K6ROT 

0 . 


PARAM 

WTMASS 

1 . 


PARAM , NOCOMPS , - 1 



PARAM 

ddrmm 

-1 


PARAM 

PRTMAXIM 

YES 



D-l 


0 


EIGRL 1 

FREQ1 2 0. .25 4096 

$ Direct Text Input for Bulk Data 
$ Elements and Element Properties for 


PSHELL 

$ Pset : 

CQUAD4 

CQUAD4 

CQUAD4 

CQUAD4 


"elprop" 

1 

2 

3 

4 


1 

will 

1 

1 

1 

1 


.00102 1 
be imported as: 
1 2 

2 3 

3 4 

4 5 


region 


elprop 

1 


M pshel 1 . 1 " 

31 

32 

33 

34 


30 

31 

32 

33 


CQUAD4 557 1 

CQUAD4 558 1 

CQUAD4 559 1 

CQUAD4 560 1 

$ Referenced Material Records 
$ Material Record : alum 
$ Description of Material : Date: 
MAT1* 1 

* A 2763 . 


576 

577 

606 

605 

577 

578 

607 

606 

578 

579 

608 

607 

579 

580 

609 

608 


7 .3 + 10 


27-Oct -00 


Time : 15:41:56 

. 3 


GRID 

GRID 

GRID 

GRID 


0 . 

0 . 

0 

. 0127 

0 . 

0 

. 0254 

0 . 

0 

.0381 

0 . 

0 


GRID 

GRID 

GRID 

GRID 

$ Loads 

SPCADD 

RLOAD1 

LSEQ 

DLOAD 


606 

607 

608 
609 

for Load Case 


.3175 

.2 54 

0 

.3302 

.254 

0 

. 3429 

.254 

0 

. 3556 

. 254 

0 

random 




$ Displacement Constraints 


of 


SPC1 


1 

B 175 
C 407 


134 

204 

436 


1 

233 
465 

$ Displacement Constraints of 
SPC1 3 134 29 

+ D 203 232 261 

+ E 435 464 493 

$ Displacement Constraints of 
SPC1 4 235 1 

$ Displacement Constraints of Load Set : 
Spci 5 235 581 THRU 

$ Pressure Loads of Load Set : pressure 

P LOAD 4 3 1 ‘ 1 • 

4 Referenced Dynamic Load Tables 


Load Set 
30 
262 
494 

Load Set 
58 
290 
522 
Load Set 
THRU 


$ Constant Load Table 


xOsimple 


59 

88 

117 

146 

291 

320 

349 

378 

523 

552 

581 


xlsimple 
87 116 

145 

174 

319 

348 

377 

406 

551 

580 

609 



yOsimple 

29 

ylsimple 

609 


THRU 


560 


TABLED 1 1 

+ F 0 . 

1. 1024. 

1 . 

ENDT 


$ 

$ RANDOM INPUT , 

, flat spectrum 

(2048 Pa 

RMS, 0- 

1024 Hz = 

$ — 

RAND PS 1 
TABRND1 101 
+TBR001 0.0 

1 1 
LINEAR LINEAR 
4096.0 1024.0 

1.0 

4096.0 

0.0 

ENDT 

101 

$ - - 

$ DAMPING TABLE, constant over 

frequency 


$ - 

TABDMP1 400 
+DP1 0.0 

+DP2 1024.0 

CRIT 

0.02 446.0 

0.0054 ENDT 

0.02 

446.00 

0.0054 


160 .2060 dB OASPL) 


SKIP 


+TBR001 


+DP1 
+ DP2 


$ Referenced Coordinate Frames 
ENDDATA b72b05e7 


D-2 


rt d o co 


Appendix E Listing of Sample Bulk Data File fixed!03_eg.bdf’ 


$ NASTRAN input file created by the MSC MSC/NASTRAN input file 
$ translator ( MSC/PATRAN Version 8 . 0 ) on October 21, 2000 a 

$ 18:06:02. 

ASSIGN OUT PUT 2 = 1 f ixedl03_eg .op2 ' , UNIT = 12 

ASSIGN OUTPUT4 = 'egvec.dat', UNIT=15, FORM= FORMATTED 

ASSIGN OUTPUT4 = 'egveH.dat', UNIT=87, FORM » UNFORMATTED 

$ Direct Text Input for File Management Section 
$ Normal Modes Analysis, Database 
SOL 103 
TIME 600 
$ 

INCLUDE 1 dm_103 - v ' 

$ 


The remaining lines are the same as those in the normal modes bulk data file (see Appendix B). 
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Appendix F Listing of Sample Bulk Data File “fixed 106_st.bdf” 


$ NASTRAN input file created by the MSC MSC/NASTRAN input 


$ translator ( MSC/PATRAN Version 8.0 ) 
$ 16:51:35 

ASSIGN OUTPUT2 
ASSIGN OUT PUT 4 
ASSIGN INPUTT4 


on October 27, 


UNIT « 12 

15 , FORM= FORMATTED 
- 17 , FORM= FORMATTED 


= ' f ixedl06_st .op2 ' 

= ■ fj.frc 1 , UNIT = 

= ' di spl . inn 1 , UNIT 

$~Direct Text Input for File Management Section 
$ Nonlinear Static Analysis, Database 
SOL 106 
TIME 600 


file 
2000 at 


$ 

INCLUDE 1 dm_10 6 . v 1 

$ 

$ Direct Text Input for Executive Control 
CEND 

SEALL = ALL 
SUPER = ALL 

TITLE - MSC/NASTRAN job created on 27 -Oct- 00 at 16:4 
ECHO = NONE 
MAXLINES = 999999999 

$ Direct Text Input for Global Case Control Data 
SUBCASE 1 

$ Subcase name : Default 
SUBTITLE=De fault 
NLP ARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1 , REAL) =ALL 

SUBCASE 2 

$ Subcase name : Default 
SUBTITLE=Def aul t 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORT1 , REAL) -ALL 
SUBCASE 3 

$ Subcase name : Default 
SUBTITLE=De fault 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORT1 , REAL) “ALL 
SUBCASE 4 

$ Subcase name : Default 
SUBTITLE=De fault 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1 , REAL) “ALL 
SUBCASE 5 

$ Subcase name : Default 
SUBTITLE=De fault 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1, REAL) “ALL 
SUBCASE 6 

$ Subcase name : Default 
SUBTITLE=De fault 
NLPARM = 1 
SPC “ 2 
LOAD = 2 

SPCFORCES ( SORT1 , REAL) “ALL 
SUBCASE 7 

$ Subcase name : Default 
SUBTITLE=De fault 
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NLP ARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORTl , REAL) =ALL 

$ Direct Text Input for this Subcase 
BEGIN BULK 


The remaining lines are the same as those in the statie analysis bulk data file (see Appendix C). 
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Appendix G Listing .f Satnp.e Bulk Data Fi.e « fl xed«>._ S t.bdr 

„ ^ hv the MSC MSC/NASTRAN input file 
$ NASTRAN input file cr ®® t ® Y 8 0 ) on October 27, 2000 at 

S translator ( MSC/PATRAN Version 8.0 ) 

$ 16:51:35. 0 p2 ' , UNIT = 12 

ASSIGN OUTPUT2 = fixedlOl, • P = lg # F 0 RM= FORMATTED 

ASSIGN OUTPUT4 = -fl^L.frc , F 0 RM= FORMATTED 

$ Nonlinear Static Analysis, Database 
SOL 101 
TIME 600 
$ 

INCLUDE , dm_101.v 

$ . £_ y pypcutive Control 

$ Direct Text Input for Executi 

CEND 

SEALL = ALL 

SUPER = ALL . . d on 27 -Oct - 00 at 16:47:55 

TITLE = MSC/NASTRAN ]ob created 

ECHO = NONE 
SUBCASE 1 

$ Subcase name : Default 

SUBTITLE=Def aul t 

nlparm - i 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1 , REAL) =ALL 

SUBCASE 2 

$ Subcase name : Default 
SUBTITLE=De fault 
NLPARM * 1 
SPC = 2 
LOAD =2 

SPCFORCES (S0RT1, REAL) = ALL 

$ Direct Text Input for this Subcase 

begin bulk 


the — g lines me ,he same as ,hose ,, ft. s,a,,c analys, * d*. me (see A*-* 
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Appendix H Listing of Sample Bulk Data File “fixeD103_eg.bdf 


$ NASTRAN input file created by the MSC MSC/NASTRAN input 
( MSC/PATRAN Version 8.0 ) on October 21 , 


$ translator 
$ 18:06:02. 

ASSIGN OUT PUT 2 
ASSIGN OUTPUT 4 
ASSIGN INPUTT4 
ASSIGN OUTPUT 4 
ASSIGN OUTPUT 4 
ASSIGN OUTPUT 4 
ASSIGN OUTPUT 4 
ASSIGN OUTPUT 4 
ASSIGN OUTPUT4 
ASSIGN OUTPUT 4 
ASSIGN OUT PUT 4 


' f ixeD103_eg.op2 1 , 
' egvec . dat ' , UNIT = 
' f ixeD . mod 1 , UNIT = 
• displ . pr2 ' , UNIT * 
1 displ -pr3 * , UNIT = 
' displ . pr4 ' , UNIT * 
' displ . inn UNIT = 
• di spL . inn ’ / UNIT a 
' egveC . dat 1 , UNIT = 
• egveH . dat 1 , UNIT j 
1 egveCa . dat ’ , UNIT 


UNIT = 12 

15 , FORM=UNFORMATTED 
65 , FORM = FORMATTED 
7 9 , FORM= FORMATTED 

80 , FORM= FORMATTED 
78 , FORM= FORMATTED 

81 , FORM= UNFORMATTED 

82 , FORM=UNFORMATTED 
11 , FORM= UNFORMATTED 
87 y FORM= UNFORMAT TED 

=> 83, FORM= FORMATTED 


$ Direct Text Input for File Management Section 
$ Normal Modes Analysis, Database 


file 
2000 at 


SOL 103 
TIME 600 


$ 

INCLUDE 1 dm_103D . v ' 

$ 


The remaining lines are the same as those in the normal modes bulk data file (see Appendix B). 
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Appendix I Listing of Sample Bulk Data File “fixeD106_st.bdf” 


$ NASTRAN input file created by the MSC MS C / NASTRAN input 
$ translator ( MSC/PATRAN Version 8 . 0 ) on October 21 , 

$ 16 : 5 1 : 3 5 . 

ASSIGN OUTPUT 2 * • fixeDl06 st.op2', UNIT = 12 

OUTPUT4 = ' f _N . f rc *7 uNIT = 15 . FORM=UNFORMATTED 
ASSIGN INPUTT4 = • displ . inn UNIT = 17 , FORM= UNFORMATTED 

$ Direct Text Input for File Management Sectron 
$ Nonlinear Static Analysis, Database 
SOL 106 
TIME 600 


file 
2000 at 


$ 

INCLUDE ' dm_106D.v' 

$ 

$ Direct Text Input for Executive Control 
CEND 

SEALL = ALL 

SUPER = ALL _ , - cc 

TITLE = MSC/NASTRAN job created on 27-Oct-OO at 16:47.55 

ECHO = NONE 
MAXLINES = 999999999 

$ Direct Text Input for Global Case Control Data 
SUBCASE 1 

$ Subcase name : Default 

SUBTITLE=Def aul t 

NLPARM - 1 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1 , REAL) =ALL 

SUBCASE 2 

$ Subcase name : Default 
SUBTITLE=De fault 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORTl , REAL) *ALL 
SUBCASE 3 

$ Subcase name : Default 
SUBTITLE= Default 

NLPARM = 1 
SPC * 2 
LOAD * 2 

SPCFORCES (SORTl , REAL) =ALL 
SUBCASE 4 

$ Subcase name : Default 

SUBTITLE= Default 

NLPARM * 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORTl , REAL) =ALL 
SUBCASE 5 

$ Subcase name : Default 
SUBTITLE=De fault 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORTl , REAL) =ALL 
SUBCASE 6 

$ Subcase name : Default 
SUBTITLE =De fault 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORTl , REAL) -ALL 
SUBCASE 7 

$ Subcase name : Default 
SUBTITLE=De fault 
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NLP ARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1 , REAL) =ALL 

$ Direct Text Input for this Subcase 
BEGIN BULK 


The remaining lines are the same as those in the static 


analysis bulk data file (see Appendix C). 
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Appendix J 


Listing of Sample Bulk Data File “fixeD101_st.bdf’ 


S NASTRAN input file created by the MSC MSC/NASTRAN input 
$ translator ( MSC/PATRAN Version 8.0 ) on October 27. 

$ 16:51:35. - 

ASSIGN OUTPUT2 = ' f ixeDlOl at.op2' f UNIT - 1 

ASSIGN OUTPUT4 = ■ £ 1 L. fret UNIT = 15 , FORM=UNFORMATTED 

ASSIGN OU , , , ttnit ”* 1? FORM= UNFORMATTED 

ASSIGN INPUTT4 = 'dispL.inn ,UNIT - 1 

S Direct Text Input for File Management Section 
$ Nonlinear Static Analysis, Database 

SOL 101 

TIME 600 


file 
2000 at 


$ 

INCLUDE 'dIn_101D.v , 


$ 

$ Direct Text Input 


for Executive Control 


CEND 

SEALL = ALL 

TITLE = MSC /NASTRAN job created on 27 -Oct -00 at 16:47:55 

ECHO = NONE 
MAXLINES = 999999999 

$ Direct Text Input for Global Case Contr 
SUBCASE 1 

$ Subcase name : Default 
SUBTITLE=De fault 


NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES (SORT1 , REAL) =ALL 

SUBCASE 2 

$ Subcase name : Default 
SUBTITLE=Def ault 
NLPARM = 1 
SPC = 2 
LOAD = 2 

SPCFORCES ( SORT1 , REAL) =ALL 

$ Direct Text Input for this Subcase 

begin bulk 


The remaining lines are the same as those in the static analysis bulk data file (see Appendix C). 
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